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Abstract 

0^ ■ For three positive integers Oj, aj, pairwise coprime, we present an algorithm that 

I find the least multiple of that is a positive linear combination of aj , ■ The average 

' running time of this algorithm is 0{1). Using this algorithm and the Chinese remainder 

theorem leads to a direct computation of the Frobenius number f(ai, 02, 03). 
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1 Introduction 

We begin with some facts and notation [2, 5]. Let ai < 02 < ^3 be relatively prime positive 
integers. It is well known that if n is sufficientely large integer the Diophantine equation 

xitti + X2a2 + xstts = n (1) 

has solutions in positives integers xi,X2, X3 and we say in this case that n is N— representable 
by ai, 02,03. Let f(ai,- ■■ , a„) be the largest integer not representable by the form (1) in 



> ■ positive integers. 

00 

o 

O 



There are representation of f(ai, 02, 03) for each i = 1, 2, 3 of the form [5]: 

.^.yij^j yij>0 (2) 

o , 

Q\ , If ai, a2, 03 are pairwise prime, we prove that in each representation above, there is only one 

term yijttj that is N— representable by ai,ak. 

We give a simple algorithm that find the least multiple of Oj that is representable by aj, a^. 
^ ■ This method is very simple (few lines of code) and similar to euclidean algorithm. It compute 

f(ai, 02, as) for large number, we test it for numbers having thousand digits and the average 
running time is 0(1) for pairwise prime triplet of numbers. 

There are other algorithms to compute f(ai,a2,a3). The Davison algorithm [3] is based on 
the LLL lattice method. The elbows algorithm is [4] is based on the integer linear program- 
ming (ILP) method. 

In section 2, we give semi-explicit representation of f{ai, aj, ak) for pairwise coprime inte- 
gers {ai,aj,ak)- We present in Section 3 our algorithm for computing the least multiple 
N— representable and prove that it terminates. 

Section 4 is devoted to randomized testing and average time needed for computation of 
f(ai, 02,03) 
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2 Representation of f(a^, a^, ak) 



For any ordering k) of the set {1, 2, 3}, let (oj, Oj, a^) three integers relatively prime. It 
is no restriction to assume that {ai,aj,ak) are pairwise coprime and none of {ai,aj,ak) is 
N— representable by the two others [2]. 

Lemma 2.1. Let 

f(ai, aj, ak) = XjUj + x^afc, Xj,Xk > (3) 

and suppose that 

Xjaj = Ujtti + Vjak, Uj > 0, vj > (4) 
Then xjaj is the smallest multiple of aj that is N— representable by ai,ak 
Proof. Suppose that f{ai,aj,ak) have another representation: 

f(ai, 02, as) = yjOj + ykttk, yj > xj, zj > (5) 

And 

yjOj = SjOi + tjOk, ^jjtj > (6) 
Combining the two expressions, we can represent f(aj, a-,-, a^) as: 

f{ai, aj, ak) = ujOj + {yj - Xj)aj + {vj + tk)ak, uj > 0, {vj + tk) > 0, vj + tk > (7) 

so f{ai,aj,ak) would be representable by 01,02,03, a contradiction to the definition of 
f{oi,aj,Ok). □ 

Lemma 2.2. For each i let f{ai,aj,ak) = VijCij + yikOk with {i,j,k) a permutation cycle of 
(1,2,3). Then either yij a j is representable by {ai,ak), or y^kOk representable by {ai,aj) 

Proof. Reresent f{ai,aj,ak) as: 

f(ai, 02, as) = ^1202 + yistts = y2iai + 7/2303 = Vsiai + 2/32^2, Uij > 0, i,j = 1, 2, 3 (8) 

Consider the couple (7/13,2/23)5 observe first that yis 7^ 2/23, otherwise 2/1202 = 2/21^1 a contra- 
diction to ai, a2 relatively prime. 
If 2/13 < 2/23 then 

2/12^2 = 2/21^1 + (2/23 - 2/13)03, 2/12 > 0, (2/23 - 2/13) > (9) 
So 2/12O2 is representable by (01,03). Alternatively the case 2/13 > 2/23 lead to 

2/21^1 = 2/12^2 + (2/13 - 2/23)03, 2/13 > 0, (2/13 - 2/23) > (10) 
And then 2/21 oi is representable by 02,03. 

We obtain the other results by considering the couples (2/12,2/32) and (2/13,2/31) n 
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Consider now integers n that is not N— representable by ai, 02, 03 but is N— representable 
by any couples {aj,ak): 



n 



Ui2a2 + Ui2a2 = U2iai + W23«3 = Usitti + ti32a2 



In the first right expression above, suppose Ui2a2 N— representable by ai, 03. Lemma 2.2, im- 
phes that the N— representable term for the third expressions must be usiCi, and consequently 
the term N— representable for the second expression «2iai. This give us a first decomposition 
of n. An other decomposition is given when we suppose ^1303 N— representable in the first 
right expression. Finally, such integers n have exactly two possible decompositions. 

To compute f (a, b, c) consider the two congruent systems: 
X = uutti moda2 ( x = ^2102 modai 

X = ^2302 mod 03 < X = uistti mod 03 

X = u^itts modai x = ^^3203 mod 02 

Solutions of the two congruent systems are given by the Chinese remainder theorem and 
f (a, b, c) is maximum of the two solutions. 



3 Least multiple of b N— representable by a, c 

Algorithm (Find least multiple representable). We are given three integers a,b,c pairwise 
prime. This algorithm find the least multiple of b that is N— representable by (a, c). 

1 . [Initialize] 

Po = -(b + c( — mod a)) 

a c 

b = Pott mod c 

2. [Sequences chains] 

Develop the integer sequences (/c,), {pi, (fj), where 



ki = 1 



c 

.Po. 



ki = l + 



Pi-2 

,pi-l 



Pi = hPo mod c, Pi = kiPi_i modpi_2, i ^ 2 



(12) 
(13) 



3. [Main loop] 

While (pi a > Vib) { 
i = i + l 

} 

return ua 



Pi ( — mod c) mod < 



i ^ 1 



Po 



(14) 
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The theoretical justification of this algorithm is as follows: 

Lemma 3.1. The integer sequence (pi) defined by equation I4 is strictly decreasing 

Proof. For i = 1, the integer ki defined by 14 is the smallest k verifying: kpo > c. In fact, 
we have from 13: pi = kipo mod c, so necessarily pi < po, otherwise {ki — 1) would be the 
least integer k satisfying kpo > c. 

For i ^ 2, the integer ki defined by 14 is the least integer k such that kpi-2 > Pi-i and with 
Pi = A;jPj„i mod Pi _2, necessarily pi < Pi-i, otherwise {ki — 1) would be the least integer k 
satisfying kpi^2 > Pi~i- 

□ 

4 Example 

Let a = 7523, b = 8231, c = 9533 



Step 


ki 


Vi 


Pi 


Pia - v,b 
c 


1 


2 


9533 


7001 


-3525 


2 


2 


7001 


4469 


-1526 


3 


3 


4469 


1937 


-1053 


4 


2 


1937 


1342 


-580 


5 


2 


1342 


747 


-107 


6 


5 


747 


152 


45 



Table 1: The least multiple of b representable by a, c is 7476 = 152c + 45a 

We test the number of step required to compute the three least multiple representable. 
We use always the scheme xb = ac — f3a, ya = r]c — ub, zc = fib — ■ya 

An implemented demonstration of this algorithm using Mathematica software is available in 
[6]. 

Compared to Davison algorithm [3] we noted that our algorithm don't need to go through 
the entire euclidean chain until reaching zero but rather the first Pi smaller than Vi. This 
seems giving substantial improvement of the running time. 
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